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We report numerical simulations on granular shear flows confined between two flat but frictional 
sidewalls. Novel regimes differing by their strain localization features are observed. They originate 
from the competition between dissipation at the sidewalls and dissipation in the bulk of the flow. 
The effective friction at sidewalls is characterized (effective friction coefficient and orientation of the 
friction force) for each regime, and its interdependence with slip and force fluctuations is pointed 
out. We propose a simple scaling law linking the slip velocity to the granular temperature in the 
main flow direction which leads naturally to another scaling law for the effective friction. 

PACS numbers: 47.57.Gc, 45.70.Mg, 83.50.Ax, 83.80.Fg 


Recent progress has been made in the theoretical de¬ 
scription of the rheology of granular materials (e.g. [l|- 
[^). However, a scientific bottleneck still prevents the 
use of those theories in real configurations: the modeling 
of the interaction of the granular material with a solid 
boundary and its implication on the rheology of the sys- 


idary ar 

tern [^, [l3- Yet, such interactions are crucial to under¬ 
stand shear banding that is a widespread phenomenon 
in slow granular flows. They are also essential for the 
comprehension of many industrial and agricultural ap¬ 
plications (silos, hoppers, chutes, mixers, and blenders) 
and geophysical phenomena, for example, the mobility of 
avalanches or the cratering due to an impact. 

The aforementioned difficulty originates from experimen¬ 
tal and numerical evidence that questions the modeling 
of a solid interface as a simple boundary condition. In¬ 
deed, the existence of cooperative effects miiii in the 
force network and in the velocity field prevents the use of 
a purely local approach, i.e., an approach in which the lo¬ 
cal stresses are simply related to the local shear rate. An 
important source of cooperative effects and thus of non¬ 
locality is the mechanical noise of the flow itself [i, [Ti¬ 
ll^. Thus, the effective friction coefficient of a granular 
material flowing on a flat but frictional interface and the 
corresponding slip velocity are partially controlled by the 
shear-induced fluctuations of the force network |16l4l8j. 
and their prediction remains a challenge. The under¬ 
standing of such interfacial phenomena is of paramount 
importance not only for granular flows but also for other 
complex fluids displaying wall slip and nonlocal effects 
such as emulsions [l9|, dense suspensions 13,[Hi, foams, 
etc. Here, by means of discrete element simulations, we 
study the effective friction of a dense granular material 
confined in a shear cell. We show that at flat frictional 
walls, even if the system is globally below the slip thresh¬ 
old, force fluctuations trigger slip , leading to a nonzero 
slip velocity and an effective wall friction scaling with a 
sliding parameter. These results shed light on the neces¬ 
sity to introduce fluctuations in theories aiming to cap¬ 
ture and predict the behavior of dense granular flows at 
the vicinity of an interface. 


Numerical simulations .— We carry out and analyze 
numerical simulations of dense granular flows in a 3D 
wall-bounded geometry. Simulations are performed by 
using the contact dynamics method [lH , as implemented 
in the LMGC90 open source framework [23|. The flow 
configuration [sketched in Fig. [Il[a)] is a rectangular 
cuboid (length lx = 20d, depth ly = lOd, and variable 
height Iz) characterized by a periodic boundary condi¬ 
tion in the main flow direction (x), two bumpy walls at 
the top and at the bottom, and two lateral flat but fric¬ 
tional walls (normal to the y direction). Gravity acts on 
the system along 2 :, and the flow is driven by the mo¬ 
tion along X of the bottom wall. The top wall cannot 
move on the x and y directions but is free to move in 
the 2 : direction, simply according to the balance between 
its weight and the force exerted by the grains. Simula¬ 
tions were performed with N = 10000 slightly polydis- 
perse spheres (uniform number distribution in the range 
O.Od — l.ld) interacting through perfectly inelastic col¬ 
lisions and Goulomb friction (/i = 0.5). The coefficient 
of restitution has nearly no influence on dense granular 
flows due to the presence of enduring contacts [2^ . Con¬ 
sequently, we chose perfectly inelastic grains to maximize 
dissipation and thus save computation time. Interactions 
of particles with the flat walls were also perfectly inelas¬ 
tic and frictional (with a coefficient of friction /ipw)- We 
performed several simulations varying the velocity of the 
bottom wall V, the weight of the upper wall M, and 
the particle-wall friction coefficient jUp^. The first two 
parameters can be made dimensionless, for example, by 
considering a particle Fronde number V = V/^/^ and 
the ratio between the mass of the top wall and the mass 
of the grains, M = MjNvn, where m is the average parti¬ 
cle mass. This flow configuration is interesting for study¬ 
ing wall friction, since slip velocity and fluctuations are 
less dependent on the wall properties than if shear acted 
on planes parallel to the wall. This does not mean that 
the wall properties do not affect velocity profiles, but that 
the problem is a little more conceptually decoupled than, 
for example, in inclined chute flows. Preliminary analysis 
indeed showed that velocity and velocity fluctuation pro- 
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FIG. 1: (color online) (a) Sketch of the wall-bounded flow 
configuration, (b) Velocity profiles at the sidewalls along z 
obtained for V = 10, M = 0.2, and different values of the 
particle-wall friction coefficient: /Xpw = 0.05, 0.1, 0.2 and 0.3. 
Three patterns of localization appear, (c) Fluctuation of the 
X particle velocities at the wall along z for the same set of 
parameters. Inset: Ratio of vertical to horizontal fluctuations 
versus z. 


files were nearly uniform in the y direction for the range 
of parameters considered in this study (V = 0.1 — 10, 
M = 0.2 — 2, /ipw = 0.05 — 0.3). The profiles at the 
wall can therefore be taken also as a reference for the 
internal behavior; shear principally acts on a plane or¬ 
thogonal to the wall. It is clear that not all variables are 
uniform in y\ for example, solid fraction 0 is influenced 
by the presence of walls and displays fluctuations around 
a decreasing value when approaching the wall [IBl. On 
the other hand, due to wall friction, stress components 
(jyz and (jyx (and their symmetric counterparts) will vary 
with y. In this work, we focus mainly on wall behavior, 
postponing a full analysis of the behavior of the system 
including also profiles in the y direction to a more de¬ 
tailed study. 

We computed average profiles along z by performing 
averages on slices with a thickness 2d in the vertical 
direction. The streamwise effective wall friction coeffi¬ 
cient at the lateral walls is estimated for each slice as 
the ratio of the average force in the flow direction x and 
the average force in the direction normal to the wall y: 
fiw = (Fx)/(Fy) [l^. This corresponds to the stress ra¬ 
tio (Jyx/cjyy at the wall. Profiles of velocity fluctuations 
in the x and 2 : direction (which are related to granular 
temperature) are also calculated as Tx = ^{vx — (t’a^))^), 
where the fluctuations are correctly 
computed with respect to the avep^e velocity value ex¬ 
trapolated at the particle center [26|. Concerning bulk 
profiles, we computed for each slice the average value 
of solid fraction (j) (not discussed here) and stress com¬ 


ponents (Txx^ o-yy^ (Jzz^ and (Jxz- The latter allow to 
define a bulk friction coefficient fixz = ^xzlv where 

P — {^XX F CFyy + (^zz)l‘^' 

Flow profiles .— An example of the velocity profiles 
obtained in this configuration is given in Fig. [IJb), for 

V = 10, M = 0.2, and different values of the particle- 
wall friction coefficient. Profiles typically display shear 
localization. For the range of parameters considered in 
this study, three regimes appear: (A) for high M and/or 
high wall friction, shear is localized at the bottom; (B) 
for low M and low wall friction, shear is localized near the 
top; and (C) for low M and intermediate wall friction, 
a central plug can form with two shear zones near the 
bumpy walls. Concerning such profiles, in the following 
we will refer to the “shear zone” as the zone where most 
of the velocity variation occurs (say, 99%) and to the 
“creep zone” as the nearly stationary remaining part. For 
regimes B and C, far from the shear band, it is better to 
refer to a “plug flow zone.” In shear zones, velocity pro¬ 
files are characterized by an exponential variation whose 
characteristic length is mainly a function of and M. 
If M ^ 1, bottom localization seems to be more proba¬ 
ble and the characteristic length decreases with //pw; it 
is reasonable to infer that a linear profile should be ob¬ 
tained for sufficiently large M and sufficiently low //pw 
Such a coexistence of different localization patterns was 
recently reported for a different flow configuration [^ . 
Figure [TJc) displays sample profiles of streamwise veloc¬ 
ity fluctuations Tx. These agree qualitatively with shear 
rate profiles, with higher fluctuations in the shear zones. 
Vertical fluctuations, as shown in the inset in Fig. [U^c), 
tend to be more important than streamwise ones in the 
creep zone in regime A. In the shear zones, there was also 
evidence of the formation of secondary (weak) convective 
rolls, similar to those observed in Refs. 00. These 
rolls appeared to be confined to the shear zone, to be 
stronger for regime C (where the shear zone was larger), 
and were characterized by particles moving upwards near 
the flat walls and descending in the center of the channel. 
The localization pattern observed, as well as the convec¬ 
tive instability, is very new and interesting and will be 
addressed in a detailed study. 

Effective wall friction .— Figures [2f a) and [2fb) dis¬ 
play profiles of the effective wall friction coefficient, for 

V = 10, M = 0.2, and different values of //pw It is 
clear that, as was seen for 2D flows [lH or in 3D con¬ 
fined gravity-driven flows [s^, Pw/Ppw decreases with 
/ipw These profiles are strongly related to velocity pro¬ 
files: in regime B, far from the shear band, grains move 
as a plug in the x direction with velocity V. There are no 
stick events, so the effective coefficient of friction in the 
plug zone is equal to /ipw In the shear zone, stick-slip 
events may emerge, and therefore < Ppw Hi- For 
the simulations where shear is localized at the bottom 
(regime A), we find a decreasing profile of y^ versus 2 : 
in the shear band and a less pronounced variation in the 
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FIG. 2: (color online) Profiles along z of (a) the effective wall 
friction coefficient fiw and of (b) fiw rescaled on the particle- 
wall friction coefficient /Xpw, for V — 10, M = 0.2, and differ¬ 
ent values of //pw = 0.05,0.1,0.2,0.3. jJLw/lJipw decreases with 
/Xpw. (c) Prohles of the bulk friction coefficient fixz — cFxzjp 
along z. The labels are the same as those used in Fig. [TJb). 


creep zone. This can be explained by the fact that, rea¬ 
sonably, stick-slip events become more and more proba¬ 
ble when we approach the creep zone M- When a central 
plug forms (regime C), a more complex profile may 
appear, displaying a flat local maximum corresponding 
to the plug zone. These results confirm what was al¬ 
ready known for 2D flows, that is, the inte rdep endence 
of effective friction, slip and fluctuations [171, ll 8 | . The ef¬ 
fective bulk friction coefficient fixz as shown in Fig. [2fc) 
qualitatively follows the behavior of velocity fluctuations 
and shear rate. An increase in wall friction contributes 
to weakening of the effective bulk friction in the creep 
zone. The effective bulk friction coefficient being gener¬ 
ally higher than the effective wall friction coefficient, in 
the creep zone they can take similar values. 

Force orientation at the wall .— The probability den¬ 
sity function of tangential force orientation depends on 
the flow regime and flow zone. For regime B, in the plug 
flow zone, tangential forces are purely oriented along x 
(counterflow), and in the shear zone only slight differ¬ 
ences exist. For regimes A and C, in the shear zone [see, 
for example, curves at 2 ; = 3d in Figs. [3l[a) and mb)], 
the most likely direction is again the x direction, but sig¬ 
nificant vertical fluctuations are present. The likelihood 
of such transversal fluctuations increases with the wall 
friction. For these two regimes, the probability of finding 
vertical force fluctuations is more important for grains 




FIG. 3: (color online) Polar representation of the probability 
density function of tangential force orientation at the wall, 
for y = 1, M = 0.2, at several heights z (= 3d, 9d, 18d, and 
36d), for (a) /Xpw =0.1 (regime G) and (b) /Xpw =0.3 (regime 
A). Statistics are performed on a 5d-wide slice centered on 
the given z value. 


close to the top stationary wall. In the creep flow zone of 
regime A [curve at z = 36d in Fig.^b)], we see that ver¬ 
tical force orientations can become more probable than 
horizontal ones. This is obviously related to the incidence 
of vertical velocity fluctuations discussed previously. The 
average tangential force is always directed opposite to 
the motion, but a very wide distribution of orientations 
is present. The probability distribution function (PDF) 
is roughly ellipsoidal. For regime C in the plug flow zone, 
due to the partial slip of the plug, the PDF is similar to 
the creep zone PDF but is bent towards the x direction, 
assuming a sort of “bunny ears” shape. Note that the re¬ 
ported behavior is similar to that obtained for confined 
granular flows on a heap, where a modification of the ori¬ 
entation of the friction force which is significant in the 
creeping zone has been reported [l^ . 

Wall friction scaling .— In order to study more in de¬ 
tail effective wall friction, we discuss now its scaling on 
system parameters. For 2D dense flows on an incline, it 
was shown that //ipw scales with a dimensionless slip 
length defined as Vs\ip/^d [TH. The origin of this scaling 
was postulated to be the presence of force and velocity 
fluctuations, yielding stick-slip events. The inverse of 
the shear rate 7 “^ was there used as an estimate of the 
time scale of force fluctuations. In the 3D flow config¬ 
uration discussed here, this scaling has some drawbacks 
and seems not to hold well in the creep zone. One reason 
for this behavior is that jd may not be a good veloc¬ 
ity scale for slip events in creeping flows, as seems to be 
the case for dense flows. Another estimate for a velocity 
scale related to force fluctuations may be the correla¬ 
tion of velocity fluctuations, e.g., granular temperature. 
The latter quantity was indeed used as a velocity scale 
in boundary conditions for flat frictional walls in the ki¬ 
netic theory of granular gases Our simulations 

show that the square root of the velocity fluctuations in 
the X direction globally scales with the shear rate 

(see Fig. SI). However, often in the creep or plug zone. 
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FIG. 4: (color online) Comparison between a velocity scale 
based on shear rate, ddvxjdz^ and one based on a:-velocity 
fluctuations, \fT^ for M = 0.2 and V = 0.1 (squares), 
M = 0.2 and V = 1 (diamonds), and M = 0.2 and y = 10 
(downward triangles), M = 2 and y = 1 (circles). The colors 
(gray scales) correspond to different grain-wall friction coef¬ 
ficients [see Fig. [Ub)]. The fluctuational velocity scale dis¬ 
plays a plateau for each simulation for low shear rates. Inset: 

increases when decreasing inertial number I. The 
labels are the same as those used in Fig. [T][b). 


when 7 decreases, the fluctuations decrease less rapidly, 
eventually displaying a plateau. This is evident from the 
inset in Fig. 01 where is shown to increase when 

decreasing the inertial number / = yd/ p (creep is 
found for / ^ 0). The reason for this plateau can be 
related to the nonlocality of the rheology: it was already 
shown that fluctuations propagate far from the shear 
Such nonlocal effects can be modeled 


zone 


UM 


by the diffusion of ffnctuating energy, as in kinetic theo¬ 


ries of dense granular flows 


, or by fluidity-based 


models [l[. Note that the plateaus are even more evi¬ 
dent for vertical velocity fluctuations (v^? iiof shown), 
which are probably influenced by vertical fluctuations of 
the position of the upper wall. 


The correspondence between the plateau of the curves 
^/T^ versus ddvxjdz and the jiyj profile in the creep 
zone suggests to test the \fT^ scaling for the slip ve¬ 
locity. This test is shown in Fig. [5l where we display the 
rescaled effective friction coefficient //ipw as a function 
of Vxl\fT^. Note that it is preferable to use the stream- 
wise velocity scale \/7k (and neither nor \/Tx + T^) 
given that also is based on the forces oriented along 
the main flow direction. Note also that Ty might be a 
better choice for more dilute flows where most grains in¬ 
teracting with sidewalls bounce on them. In the present 
study, flows are dense, and grains mostly slide on side- 
walls. We can see that the scaling performs globally well 
on several orders of magnitude. There are indeed some 
deviations from the master curve, in particular, some 
simulations in regime C (orange squares, diamonds, and 


FIG. 5: (color online) Rescaled effective wall friction coeffi¬ 
cient versus the dimensionless slip parameter Vx/\/T^- The 
inset shows that a power law relationship holds between 
VxlV^- The labels are the same as those 

used in Fig. 01 


triangles in Fig. 0 which show slightly higher values of 
Mw/Mpw than the other simulations. This could be re¬ 
lated to the presence of convective structures discussed 
previously. This kind of structure is known to affect ef¬ 
fective friction 28|, [ 2 ^ and may be the cause of the slight 
discrepancy between these simulations and the rest of the 
set. The functional form of the scaling appears to be the 
same as in Ref. [HI (with the exception of the different 
velocity scale): 


/ipw 


A + 


B • 


For the present case, we estimated ^4 = 1.2 and B = 1.7 
by curve fitting. It should be pointed out that this scaling 
is similar to the boundary conditions of kinetic theories 
of granular flows [Mi and used in the extension of the 
kinetic theory proposed by Jenkins and Berzi 35, Is^ . 
Moreover, a connection also exists with the nonlocal the¬ 
ories recently developed EHi, since the granular fluid¬ 
ity, i.e., the ratio of the pressure to the viscosity, can 
be shown to scale with the square root of the granular 
temperature. 

Conclusions .— Results presented in this work con¬ 
firm for 3D flows the conceptual framework which was 
suggested in Ref. 113 ; at fiat frictional walls, force fluc¬ 
tuations trigger slip events even if the system is globally 
below the slip threshold. These stick-slip events produce 
(i) a nonzero average slip velocity and (ii) a variable ef¬ 
fective wall friction coefficient which scales on a dimen¬ 
sionless slip parameter. It is interesting to see that this 
picture of wall slip (fluctuations leading to a finite proba¬ 
bility of yielding) is conceptually the same as that at the 
heart of recent findings on granular creep flows Hi Hi- 
The scaling of the slip velocity on the rms velocity fiuctu- 
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ations supports also the idea that fluctuations are a key 
ingredient that should be included in any model aim¬ 
ing to describe dense granular flows at the vicinity of an 
interface. The scaling law can be used as a boundary 
condition in a nonlocal theory for dense granular flow 
implying granular temperature 35,33- In order to solve 


such a model, an additional constraint is needed: a pos¬ 
sible choice is to equate the flux of fluctuating energy at 
the boundary to the dissipation, as in kinetic theories. 
An expression for the energy dissipation based on fric¬ 
tional sliding should be studied. In Ref. 0, a scaling 
based on the shear rate was shown to be effective for 2D 
inclined chute flows. In that case, given that the shear 
acted in planes parallel to the wall, the scaling law could 
be interpreted as a Navier partial slip law. The present 
contribution shows that such a partial slip boundary con¬ 
dition involving the velocity gradient normal to the wall 
cannot be universal. Fluctuations trigger slip events and 
therefore affect effective friction, whether shear planes 
are parallel or not to the wall. Thus, if a scaling on shear 
rate was valid in that case, it was probably because the 
shear rate was a good measure of the scale of fluctuations 
for the considered data. 

Understanding the behavior of granular materials at a 
flat but frictional interface is a preeminent scientific chal¬ 
lenge. In this work, we demonstrated the relevant role 
of fluctuations and pushed forward our understanding of 
effective wall friction and wall slip. In the future, it will 
be interesting to test the framework relating fluctuations, 
wall friction and wall slip in other contexts, the problem 
of interface rheology being common to other divided me¬ 
dia, such as foams, clays, or gels. The effect of grain 
shape should also be explored, since it may influence the 
wall-slip velocity and thus the reported scalings. 

The numerical simulations were carried out at the 
CCIPL (Centre de Calcul Intensif des Pays de la Loire) 
under the project MTEEGD. 
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